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Abstract 

Bose-Einstein condensation (BEC) of Feshbach molecules in a homogeneous Bose gas is studied 
at finite temperatures in a single-channel mean-field approach where the Hartree-Fock energy and 
pairing gap are determined self-consistently. In the molecular-BEC state, the atomic excitation 
is gapped and the molecular excitation is gapless. The binding energy of Feshbach molecules is 
shifted from the vacuum value due to many-body effect. When the scattering length a s of atoms is 
negative, the system is subject to mechanical collapse due to negative compressibility. The system 
is stable in most regions with positive scattering lengths. However at low temperatures near the 
resonance, the molecular-BEC state vanishes, and the coherent mixture of atomic and molecular 
BEC is subject to mechanical collapse. 



* Electronic address: yinlan@pk u.edu.cn 



1 



I. INTRODUCTION 



Many interesting phenomena have been observed in Bose gases with Feshbach resonances, 
normous particle-loss rate due to three-body recombination was found near the resonance 
lfl . Oscillations between atoms and diatomic Feshbach molecules were generated by a sudden 

n 

change in magnetic field [2[. Large number of Feshbach molecules were produced from Bose 

ating the 
4j. Even 



gases either by tuning the magnetic field through the resonance [3| or by oscil 



magnetic field at a frequency corresponding to the molecu 



tetramer molecules were created from Feshbach molecules {5]. However BEC of Feshbach 
molecules in Bose gases has not been achieved in experiments so far and properties of this 
state are waiting to be explored. 



ar binding energy 



Properties of atomic-BEC state have been extensively studied theoretically js]. Phase 



transition between atomic and molecular BEC was proposed near the resonance j7-ll0|. 



However these phases suffer mechanical collapse in regions with negative scattering length 



a s < A recent work 12] shows that the molecular-BEC state of a Bose gas with a 
wide Feshbach resonance at zero temperature exists only when the atom density n satisfies 
na^ < 0.0164 for positive scattering length a s > 0. 

In this paper, the molecular-BEC state of a homogeneous Bose gas with a wide Feshbach 
resonance is studied at finite temperatures. In the following, we first describe the molecular- 
BEC state in a mean-field approach where the Hatree-Fock energy and pairing gap are 
determined self-consistently. Long-wavelength excitations and molecular binding energy are 
obtained. Then mechanically stability of the molecular-BEC state is examined across the 
resonance, and the mean-field phase diagram is obtained. Coherent mixture of atomic and 
molecular BEC is also studied. Discussion and conclusion are given in the end. 

II. MEAN-FIELD THEORY OF THE MOLECULAR-BEC STATE 
A. Mean-field approach 

In alkali-atom gases with Feshbach resonances, scattering states in the open channel 
and bound states in the closed channel are coupled together. Near a wide resonance, the 
effective range of the interaction is very small, and most atoms are in the open channel 

nn 

[131 . Il4| . A uniform Bose gas with a wide Feshbach resonance can be effectively described 



by a single-channel model, 

fi 2 

where tp is the field operator of Bose atoms in the open channel, the coupling constant is 
given by g — 4irh 2 a s /m, and /i is the chemical potential of atoms. 

In the molecular-BEC state, atoms are paired into molecules and the off-diagonal long 
range order (ODLRO) appears, A = g(ijjip). Without losing any generality, A > is 
assumed in the following. In the mean field, in addition to the pairing energy, both Hartree 
and Fock energies should be taken into account, and the mean-field Hamiltonian density is 
given by 

U p = H° p - |Vv 2 ?A + (2gn - + y + W), (2) 

where the constant is given by = — (gn 2 + A 2 /2g). In this mean-field approximation, 
the pairing gap A and Hartree- Fock energy are determined self-consistently, but the fluctua- 
tion and three-body effects are ignored. The mean-field theory provides a description about 
single-particle excitations, but it is incapable to describe collective excitations or few-body 
properties such as Effimov effect. Nonetheless, the mean-field approximation can serve as a 
starting point of the theoretical description for the molecular-BEC state. 

The mean-field Hamiltonian can be diagonalized by Bogoliubov transformation, 

H p = £ + Y^E k a{a^ (3) 

k 

where = u k ipk + ^fc^l^ is the annihilation operator of atomic quasi-particles, E k = 
a/ (e^ — fi') 2 — A 2 is the quasi-particle energy, and the transformation coefficients are given 
by u\ = v\ + 1 = [1 + (e fc — fi')/E k ] [2 with e k = h 2 k 2 /2m and fjf — fi — 2gn. The energy 
constant £q is given by 

k 

where 1//3 = k^T and V is the volume. 

The order parameter A can be determined self-consistently, 

A = ~^2u k v k [l + 2f k ], 

k 



E, 



A ., ff SE k , 
coth( 



2E, 



- gn 2 v, 
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where f k is the Bose distribution function of quasi-particles, f k 
self-consistency equation can be written explicitly as 



l/[exp(PE k ) - 1]. This 



■in 



Aixh 2 a* 



d 3 k 



1 (PE k 
coth 



1 

2T k 



(4) 



2E k V 2 

where the term — l/2e k in the integrand on the right-hand side is a counter term due to 
vacuum renormalization. The chemical potential \x and the order parameter A can be 
solved from Eq. fl4]) and the following equation for the density n, 

^J2H + (ui + vt) fk ] 

'(3E k \ 1 



n 



k 

d 3 k 
(2^)3 



2E h 



coth 



(5) 



B. Excitations in the long- wavelength limit 



There is a gap in the atomic excitation energy given by 



E = vV 2 - A 2 . 



(6) 



If the molecular-BEC state is stable, one necessary condition is that the gap E Q is real, or 
—jj > A. When — // = A, the atomic excitation becomes gapless, Eq = 0, which marks the 
lower boundary of the molecular-BEC phase. The upper phase boundary is located at the 
superfluid transition temperature T p where A = and Eq = —fjf. 

In addition to atomic excitations, there are collective excitations of Feshbach molecules. 
Although the molecular excitation energy cannot be obtained in the mean-field approxima- 
tion, it can be extracted from poles of two-particle correlation functions given by 



Ur[b a (r,T)bUv',r')]), 



(7) 



where b\ = ip 2 , b<i = b\, 63 = Qtbty, < r < (3h, and T is the time-ordering operator. As 



in the zero-temperature case 
phase approximation (RPA), 



12| , the correlation function can be calculated in the random- 



X (k,u;) = [l-^ )(k,u;)]-V 0) (k, 



(8) 



where ;\;( )(k, to) is the correlation function calculated in the mean-field approximation. In 
RPA, the dispersion of Feshbach molecules satisfies the following equation, 



det|I-0x (o) (k,")| =0, 



(9) 
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FIG. 1: Binding energy of Feshbach molecules in molecular-BEC state. The solid and dashed lines 
are binding energies at T = T p and T = OK. The dotted line is the binding energy in vacuum, 
Ebo — h 2 /(ma 2 ). The inset shows the temperature dependence of the binding energy for different 
n l / z a s . 

where I is the identity matrix. The molecular excitation is gapless in the long- wavelength 
limit, k = and u = 0, following Eq. (J4]). At small k and u>, to the leading order of k, the 
molecular excitation frequency is linearly dispersed. The detail of molecular excitations will 
be studied in our future work. 




C. Binding energy 

The binding energy Ef, of a Feshbach molecule in the molecular-BEC state is denned as 
the energy difference between two atomic and one molecular excitations at k = 0. Since 
the molecular excitation energy is gapless in the molecular-BEC state, the binding energy 
of Feshbach molecules is twice the energy gap of atomic excitations, Ef, = 2Eq. Due to 
many-body effect, the binding energy E b is different from its vacuum value E^ = h 2 /{ma 2 s ). 
At the superfluid transition temperature T p , the binding energy is given by E b (T p ) = —2//. 
In the limit of weakly-attractive interaction, n l ^a s — > we obtain 

E b (T p ) = lQnC 2/3 (3/2)ni\a s \ 2 k B T a , 

where ((x) is the Riemann-zeta function and T a is the ideal BEC temperature, T a = 
(2irh 2 /mkB){n/((3/2)] 2 ^ 3 ; in the limit of weakly repulsion, n 1//3 a s — > + , the binding en- 

5 



ergy is same as in vacuum, E b (T p ) = E b0 = ft 2 /(ma 2 ). 

At finite temperatures, molecules not only exist for a s > 0, but also appear as loosely- 
bound pairs for a s < in the mean-field theory, whereas in vacuum Feshbach molecules 
exist only for a s > 0. However at zero temperature, the molecular-BEC state exist only 
when l/(n 1 / 3 a s ) > 3.94 [jjj]. In Fig. [H the molecular binding energy in the mean-field 
approximation is plotted at both T = T p and T = 0, in comparison with its vacuum 
value. The temperature dependence of the binding energy for different scattering lengths 
is shown in the inset. For l/(n 1//3 a s ) < 3.94, the binding energy decreases monotonically 
with temperature, all the way to zero at the lower boundary of the molecular-BEC phase; 
for l/(n 1//3 a s ) > 3.94, the binding energy decreases monotonically to a finite value at zero 
temperature. In the dilute limit, ■n, 1 / 3 a s — > + , the temperature dependence of the binding 
energy is weak, and the binding energy is approximately given by vacuum value E b Q. 



D. Mechanical stability 



A mechanically stable Bose gas has positive compressibility, indicated by dfi/dn > 0. In 
the normal state, a Bose gas with a s < collapses when dfi/dn = 0. In the Hartree-Fock 
approximation, this condition can be rewritten as 

k B T c 1 f d 3 k 



csch^ 



2 ( e k — ^ 



(10) 



g 2 J (2tt) 3 ~ \2k B T c 
The colla pse temperature Tq is always higher than the molecular BEC transition temper 



ature T n 



15 



5Hl7j. In the limit of weakly-attractive interaction, n 1//3 a s — > , the collapse 
temperature Tq approaches to the ideal atomic BEC temperature T a , 

167T 



1 + 



3C 4/3 (| 



(11) 



In contrast, a normal Bose gas with a s > is always mechanically stable. 

We examine the compressibility of the molecular-BEC state numerically following Eq. 
fll]) and In Fig. (2J the chemical potential \i is plotted as a function of density n at given 
temperatures and scattering lengths. For a s < 0, the chemical potential fi is monotonically 



decreasing with the increase of n, du/dn < 0, showing that the system is mechanically 
unstable as found in previous studies Q- For a s > 0, the compressibility is always 
positive, dn/dn > 0, and the molecular-BEC state is mechanically stable. Therefore, in a 
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FIG. 2: Chemical potential /i versus l/n 1//3 a s . The solid line is the chemical potential in the 
molecular-BEC (MBEC) state and the dotted line is the chemical potential in the normal state, 
(a) with attractive interaction, a s < 0; (b) with repulsive interaction, a s > 0. The lines with 
IflslMdB =0.5, 0.7, 1.0, and 1.3, are plotted, where AdB = \j 'Zixfi 2 / mk-gT is the thermal de Broglie 
wavelength. 



uniform Bose gas, the molecular-BEC state can only exist with positive scattering lengths, 
a s > 0. In a trap, the molecular-BEC state may be stabilized by the finite-size effect, similar 
to trapped BEC with attractive interactions |l8| . 



E. Mean- field phase diagram 

The upper boundary of the molecular-BEC phase is determined by the disappearance 
of order parameter A = 0, corresponding to the second-order phase transition between the 
normal and molecular-BEC phases. At the transition temperature T p , Eq. (j3J) is the same 



7 



as 



the Thouless criterion which marks the pairing instability of the normal phase, 



1 

9 



d 3 k 



■ coth 



1 

2Tk 



(12) 

.2(e fe -//') \2k B T p J 2e k \ K J 

The density equation (jSJ) at T p becomes the Hartree-Fock self-consistency condition of the 
normal Bose gas 

n ~ J {2nfe^-^)/ k ^ - 1 ^ ' 

In the limit of the weakly- attractive interaction, n}^a s — > 0~, the transition temperature 
T p is close to the ideal atomic BEC temperature T a = (2irh 2 /mkB)[n/((3/2)] 2 / 3 , 

8tt 



T 



3C 4/3 (j 



-n l '*\a s 



(14) 



as found in a previous study [17| . 

In the opposite limit, n 1 / 3 a s — > + , molecules are tightly bound with the binding energy 
approximately given by E^. From Eq. (fT2|) . the chemical potential is approximately given 
by \i = —Ebo/2. The solution of density equation ( 1131) yields 



k B T p = E bo ln 1 (E bQ /k B T a )/3, 



(15) 



which is essentially the same as the mean-field result about the BEC limit of a Fermi gas 
with a BEC-BCS crossover. As pointed out in earlier works 19], the divergent mean- field 
transition temperature in this limit in fact corresponds to the temperature of molecule dis- 
sociation, not molecular condensation. The transition temperature T p can be renormalized 
by pairing fluctuations [20 1. 

The lower boundary of the molecular-BEC phase is determined by A = — //, where the 
excitation energy becomes gapless, E^ = \/ €k{ek — 2/iQ. At zero temperature, this boundary 



is located at l/(n 1//3 a s ) = 4(3/7r 



,1/3 



3.94 



121 ] . Disappearance of gap in excitation energies 



15 



201 ] . However, as discussed later 



was interpreted as appearance of atomic condensation 
in section III, a coherent mixture with atoms and molecules is mechanically unstable. 

The order parameter A in the molecular-BEC phase evolves smoothly between two phase 
boundaries. Based on the above results, the mean-field phase diagram is plotted in Fig. |3j 
The molecular-BEC phase is labeled as P(S) for a s > and P(C) for a s < 0, where symbols 
(S) and (C) represent mechanical stability and mechanical collapse. 
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FIG. 3: Mean-field phase diagram of a pairing Bose gas with a wide Feshbach resonance. P(S) 
denotes the stable molecular-BEC phase; P(C) denotes the collapsing molecular-BEC phase; N(S) 
denotes the stable normal phase; N(C) denotes the collapsing normal phase; M(C) denotes the 
collapsing coherent mixture of atoms and molecules; X denotes no solution according to the mean- 
field theory. The Dashed line is the mean-field collapse temperature of the normal phase with 
attractive interactions. 



Recently similar phase diagrams were obtained in Ref. & within the approach pioneered 
by Nozieres and Schmitt-Rink (NSR) [211 ] . In the NSR approach, pairing fluctuations are 
treated by the t-matrix approximation, and the total particle density includes not only the 
mean-field density, but also the density of thermal molecules. Although there are some 
discrepancies, most features of our mean-field phase diagram Fig. [3] qualitatively agree 



with those in Ref. 20[. The condensation temperature of the MBEC phase is significantly 
reduced in the NSR approach. In the limit of weakly-repulsive interaction, it can recover 
the condensation temperature of an ideal molecular Bose gas, T p = T a /2 5//3 . Both NSR and 
mean-field phase diagrams show that there is a transition from the molecular BEC to the 
mixture state. As shown in the next section, the mixture state is mechanically unstable in 



the mean-field approximation. In Ref. 20j, the normal state was found to be stable in the 
simple NSR approach, but unstable close to the resonance in the extended NSR approach by 
taking into account interaction between molecules . In the mean-field approach, the normal 
state is unstable near the resonance only on the attractive-interaction side. 
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III. MECHANICAL COLLAPSE OF COHERENT ATOM-MOLECULE MIX- 
TURE 



At the lower boundary of the molecular-BEC phase, the atomic excitation energy van- 
ishes, Eq = 0, suggesting the appearance of the coherent atom-molecule mixture. In the 
mixture, since both molecules and atoms exhibit ODLRO, two order parameters can be intro- 
duced 12j, ipo = ("0) describing atomic condensate and A = g{5ip5ijj) describing molecular 



condensate, where 5ip = ip — ipQ. These two order parameters are in principle independent, 
whereas in the atomic-BEC state A is a function of ipQ. In the mean-field approximation, 
the Hamiltonian density for this mixture state is given by 



Km = H° m - |-^tV 2 ^ + (2gn - fi)5^ 
2m 



(16) 



where = gn\j2 — fin — g5n 2 — \A\ 2 /2g, n = \ipo\ 2 is the atomic-condensate density, 
5n = (Sip'Sip), and n = no + Sn. The order parameter ipo should minimize the mean-field 
thermodynamic potential, which leads to the saddle point condition 



[i = g(n + 25n) + A 



(17) 



To satisfy the saddle point condition, the product A?p 2 must be real. For simplicity, we 
choose A > 0. 

As in the molecular BEC case, the mean- field Hamiltonian of the mixture can be diago- 
nalized by Bogoliubov transformation, where the field operator of atomic quasi-particles is 



given by a k = u k i(j k +v k i}jl 
= V ( e fc - A*') 



k> U k 



[l + (e k — /j/)/E k ]/2, and the quasi-particle energy 



is given by E k 
the self-consistency equation 



(gipo + A) 2 . The order parameter A can be determined by 



i.e. 



rn 



+ 



V', 



AttK a s A + g% 



d 3 k 
(2tt) 3 



2 Ei. 



coth 



(3E k 



1 

2T k 



(18) 



where similar to Eq. ((4]) the last term in the integrand is the counter term. The equation 
for the total density is given by 



n 



n + 



d 3 k 
(2tt)3 



2E h 



coth 



($E h 



(19) 
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FIG. 4: Chemical potential [i of the coherent mixture versus 1/n 1 / 3 ^ with a s /A dB =0.2, 0.3, 0.5 



and 1.0 respectively, where AdB = y2irh jvnk^T. 



From the saddle point condition Eq. ( |T71) . the ratio ipo/ipo is either 1 or -1, and the 
quasi-particle excitation energy can be rewritten as 

E k = ^{e k + 2gn ){e k - 2A^/^ )- 

If the excitation energy E k is real, two conditions, gn > and ipo/ipo = —1, must be 
true, which cannot be satisfied when the scattering length a s is negative, indicating that the 
mixture cannot exist with a s < 0. When ip = 0, there is a transition from the coherent 
mixture phase to the molecular-BEC phase. At zero temperature, the transition point 
between the two phases locates at l/(n 1//3 a s ) ~ 3.94. 

For fixed total density, the chemical potential and order parameters can be solved from Eq. 
(IT71 dU IT9]) . The solution shows a reentrant behavior near the transition to the molecular- 
BEC phase, as shown in Fig. HI The reentrance range in temperature is broadened close to 
the resonance. This reentrance behavior is likely to be an unphysical result in the mean-field 
approximation near the transition, similar to that in the Popov approximation of a Bose 
gas near the atomic BEC temperature 22| . Away from this reentrance, as shown in Fig. HI 
the chemical potential \x is a monotonically decreasing function of density n for fixed T and 
a s , dfi/dn < 0, meaning that compressibility is negative. Therefore the coherent mixture is 
subject to mechanical collapse. 
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IV. DISCUSSION AND CONCLUSION 



bck energy are 



20], the pairing 



In the mean-field approximation, although the pairing gap and Hartree- 
computed self-consistently, fluctuations are ignored. In the NSR approach 
fluctuation is treated, which significantly reduced the condensation temperature. However, 
the Hartree-Fock energy is ignored in the NSR approach, affecting the condition of me- 
chanical stability. In future many-body theories, the self-energy should be determined self- 
consistently with fluctuations properly treated. Experimentally, three-body recombination 
causes enormous particle-loss near the resonance. In a more accurate microscopic theory, 
both many-body and few-body effects should be considered. Nonetheless, the mean-field 
theory provides a basic picture about the strongly interacting Bose gas, with most features 
of the phase diagram in qualitative agreement with the NSR approach 20]. 



In conclusion, the molecular-BEC state of a homogeneous Bose gas with a wide Feshbach 
resonance is studied at finite temperatures. In the long- wavelength limit, the atomic excita- 
tion is gapped and molecular excitation is gapless. The molecular binding energy is changed 
from its vacuum value due to many-body effect and can even vanish when the scattering 
length a s is finite. When a s < 0, the compressibility of the molecular BEC state is negative 
and the system is subject to mechanical collapse. The possibility of coherent atom-molecule 
mixture is also explored, but the compressibility of the mixture is negative and the mixture 
is subject to mechanical collapse. Based on these results, a mean-field phase diagram at fi- 
nite temperatures is obtained. This work is supported by NSFC under Grant No. 10674007 
and No. 10974004, and by Chinese MOST under grant number 2006CB921402. 
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